function [nzij_pred, nzij_current, nzij_fwrd] = dynamic_g1_nz()
% Returns the coordinates of non-zero elements in the Jacobian, in column-major order, for each lead/lag (only for endogenous)
  nzij_pred = zeros(6, 2, 'int32');
  nzij_pred(1,1)=12; nzij_pred(1,2)=2;
  nzij_pred(2,1)=13; nzij_pred(2,2)=3;
  nzij_pred(3,1)=11; nzij_pred(3,2)=10;
  nzij_pred(4,1)=19; nzij_pred(4,2)=16;
  nzij_pred(5,1)=17; nzij_pred(5,2)=17;
  nzij_pred(6,1)=18; nzij_pred(6,2)=18;
  nzij_current = zeros(52, 2, 'int32');
  nzij_current(1,1)=16; nzij_current(1,2)=1;
  nzij_current(2,1)=12; nzij_current(2,2)=2;
  nzij_current(3,1)=14; nzij_current(3,2)=2;
  nzij_current(4,1)=16; nzij_current(4,2)=2;
  nzij_current(5,1)=13; nzij_current(5,2)=3;
  nzij_current(6,1)=15; nzij_current(6,2)=3;
  nzij_current(7,1)=16; nzij_current(7,2)=3;
  nzij_current(8,1)=1; nzij_current(8,2)=4;
  nzij_current(9,1)=10; nzij_current(9,2)=4;
  nzij_current(10,1)=14; nzij_current(10,2)=4;
  nzij_current(11,1)=2; nzij_current(11,2)=5;
  nzij_current(12,1)=10; nzij_current(12,2)=5;
  nzij_current(13,1)=15; nzij_current(13,2)=5;
  nzij_current(14,1)=1; nzij_current(14,2)=6;
  nzij_current(15,1)=3; nzij_current(15,2)=6;
  nzij_current(16,1)=9; nzij_current(16,2)=6;
  nzij_current(17,1)=2; nzij_current(17,2)=7;
  nzij_current(18,1)=4; nzij_current(18,2)=7;
  nzij_current(19,1)=9; nzij_current(19,2)=7;
  nzij_current(20,1)=3; nzij_current(20,2)=8;
  nzij_current(21,1)=7; nzij_current(21,2)=8;
  nzij_current(22,1)=8; nzij_current(22,2)=8;
  nzij_current(23,1)=11; nzij_current(23,2)=8;
  nzij_current(24,1)=18; nzij_current(24,2)=8;
  nzij_current(25,1)=4; nzij_current(25,2)=9;
  nzij_current(26,1)=7; nzij_current(26,2)=9;
  nzij_current(27,1)=8; nzij_current(27,2)=9;
  nzij_current(28,1)=11; nzij_current(28,2)=9;
  nzij_current(29,1)=17; nzij_current(29,2)=9;
  nzij_current(30,1)=1; nzij_current(30,2)=10;
  nzij_current(31,1)=2; nzij_current(31,2)=10;
  nzij_current(32,1)=3; nzij_current(32,2)=10;
  nzij_current(33,1)=4; nzij_current(33,2)=10;
  nzij_current(34,1)=11; nzij_current(34,2)=10;
  nzij_current(35,1)=9; nzij_current(35,2)=11;
  nzij_current(36,1)=10; nzij_current(36,2)=12;
  nzij_current(37,1)=20; nzij_current(37,2)=12;
  nzij_current(38,1)=8; nzij_current(38,2)=13;
  nzij_current(39,1)=5; nzij_current(39,2)=14;
  nzij_current(40,1)=6; nzij_current(40,2)=14;
  nzij_current(41,1)=1; nzij_current(41,2)=15;
  nzij_current(42,1)=2; nzij_current(42,2)=15;
  nzij_current(43,1)=3; nzij_current(43,2)=15;
  nzij_current(44,1)=4; nzij_current(44,2)=15;
  nzij_current(45,1)=6; nzij_current(45,2)=15;
  nzij_current(46,1)=19; nzij_current(46,2)=16;
  nzij_current(47,1)=17; nzij_current(47,2)=17;
  nzij_current(48,1)=18; nzij_current(48,2)=18;
  nzij_current(49,1)=20; nzij_current(49,2)=19;
  nzij_current(50,1)=5; nzij_current(50,2)=20;
  nzij_current(51,1)=7; nzij_current(51,2)=20;
  nzij_current(52,1)=19; nzij_current(52,2)=20;
  nzij_fwrd = zeros(4, 2, 'int32');
  nzij_fwrd(1,1)=3; nzij_fwrd(1,2)=8;
  nzij_fwrd(2,1)=4; nzij_fwrd(2,2)=9;
  nzij_fwrd(3,1)=6; nzij_fwrd(3,2)=15;
  nzij_fwrd(4,1)=6; nzij_fwrd(4,2)=20;
end
